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INTRODUCTION 


The prospect of extensive space operations in the era of the Space Transporta- 
tion System and possibly a permanent manned space laboratory accentuates the need for 
understanding the interaction of energetic heavy ions with extended matter (ref* 1). 
In addition to the space program, high-altitude aircraft operations (ref* 2) and 
radiotherapy (ref. 3) are areas requiring better understanding of heavy ion beam 
transport* In the present report, previous calculations (refs. 4 and 5) are extended 
to arrive at a more complete model for heavy ion beam transport. 

The solution of the heavy ion transport equations is considered in the straight 
ahead approximation. The equations are solved by methods of characteristics using a 
perturbation expansion. After the transport coefficients are reviewed, results for 
neon beams in water are presented along with conclusions based on these results. 


HEAVY ION TRANSPORT 

In moving through extended matter, heavy ions lose energy through interaction 
with atomic electrons along their trajectories. On occasion, they interact violently 
with nuclei of the matter producing ion fragments moving in the forward direction and 
low energy fragments of the struck target nucleus. The transport equations for the 
short range target fragments can be solved in closed form in terms of collision den- 
sity (ref. 4). Hence, the projectile fragment transport is the interesting unsolved 
problem. In previous work, the projectile ion fragments were treated as if all went 
straight forward (ref. 5). 

With the straight ahead approximation and the target secondary fragments 
neglected (refs. 4 and 5), the transport equation may be written as 



where <£j(x,E) is the flux of ions of type j with atomic mass Aj at x moving 
along the x-axis at energy E in units of MeV/amu, Cj(E) is the corresponding 
macroscopic nuclear absorption cross section, Sj (E) is the change in E per unit 
distance, and m- k (E) is the multiplicity of ion j produced in collision by 
ion k. (A listof symbols appears after the references in this report.) The range 
of the ion is given as 


Rj(E) 



dE 1 

S j ( E* ) 


( 2 ) 


The solution to equation (1) is to be found subject to boundary specification at 
x = 0 and arbitrary E as 


$j(0,E) = F j ( E) 

Usually Fj(E) called the incident beam spectrum. 

It follows from Be the 1 s theory (ref, 6) that 


A 2 . 

Sj (E) = £ - 3 -- Sp(E) 

A j Z p 


for which 



R. (E) 
D 



R p (E) 


( 3 ) 


(4) 


(5a) 


The subscript p refers to proton. Equation (5a) is quite accurate at high energy 
and only approximately true at low energy because of electron capture by the ion 
which effectively reduces its charge (ref. 7), higher order Born corrections to 
Bethe's theory (ref. 8), and nuclear stopping at the lowest energies (ref. 9). 

Herein, the parameter Vj is defined as 

v. R . (E) = v, R, (E) (5b) 

3 3 k k 


so that 


v. = Z. 2 /A. (6) 

3 3 3 


Equations (5) and (6) are used in the subsequent development and the energy variation 
in v. is neglected. The inverse function of Rj(E) is defined as 


E = Rj 1 [Rj (E) ] 


(7) 
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and plays a fundamental role subsequently. For the purpose of solving equation (1), 
define the coordinate transformation (refs. 4 and 10), 




. = x - R. (E) 

D 


= X + R.(E) 


and new functions 


Xj ( vV E E j (E) V X ' E) 


Xk ( ^ j 9 ^ j ) “ ^ ^ ^ 


where 

^ + ^ + \ 

J D 

for which equation (1) becomes 



( 8 ) 


(9) 


( 10 ) 



m jk°k 






'V 


(ID 


where the Oj are assumed to be energy independent. Solving equation (11) by using 
line integration with an integrating factor. 


wv 


= exp 


J + V 


results in 


( 12 ) 


VVV -eKp[-l<ys. + n 3 )] x^'VV 

+ dT1 ' 


(13) 
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where 


\ + v i \ _ v i ^ 

n' — ti 1 + — r 

\ 2v k 2v k S 


v k " v -i \ + v i 

= A 1 n ’ + 3. 




( 14 ) 


and the boundary condition (eq. (3)) is written as 




Consider a Neumann series for equation (13) for which the first term is 


x j 0) (i y *4 ) = exp [- i a j (n j + ^)] Sj [r- 1 (^) ] f.o. 1 ^.)] 


(15) 


and the second term is 


r^j v. 

( vy "T J . V”' • ”j > ] Z m jk\ TT exp [- T t 'k < ' 1 k + ik’l 

"5. 1 - k 


* SW'tk*] p k[ R k 1( 5k>] dT i' 


(16) 


An expression for xj^(Tj.,£j) is derived once equation (16) is reduced and higher 
order terms can be found by continued iteration of equation (13). These expressions 
(eqs. (15) and (16)) are now simplified for a monoenergetic beam of type M ions. 

The boundary condition is now taken as 


Fj(E) = 6 jM 6(E - E 0 ) (17) 

where 6j M is the Kronecker delta, 6( ) is the Dirac delta, and E Q is the inci- 
dent beam energy. Thus, 
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X)<-t)-5j> ■ 5 j[^’ ( 5j>] 6 j« - « jM «[Ej-Rj<V) 


(18) 



for which 




( 0 ) 


becomes 


*S 0) <VV - 6 j« ^p[- l ♦tj)] *[5 j - B j <B o > ] 


and ^ ^ becomes 


(1) 


s' ( w 


r* j v. 

1 I m j« “h ^ exp [j CT j <t >' - y - 1 * 4 a>] 

F . M 

D 


* ^m-VVI d i' 


The contribution to the integral (eq. (20)) occurs at 


2v m V m + V n 

' ITT R « (E °’ - ^ ^ 


provided that n' lies on the interval -£j < t] 1 < rjj, so that 


Xj 1)(i yV = exp f~ 1 0 M^j + T » ,) - i a j (T1 j ■ t > ,) 3 


M D 


The simplified form in equation (2 2) may now be used to calculate the next 
of equation (13): 


Vj /*T)j 

j 2> <TVj ,Cj ) m jk<V\M q M v M - V), jf.^ ex P[-T °M<^ + ^ 

jk j 


“ 7 °k (l Tk " ^ " 7 a j ( ’Tj " T )" > ] dT l" 


(19) 


(20) 


( 21 ) 


( 22 ) 


iteration 


(23) 
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where 


v. + v. v. - v. 

.1 k i k i „ 

t - ~ 2 ^~ ' + -T^~ 

k k 


v, - v. v, + v. 

JS 1 


■> 


( 24 ) 




2 \ l i 




and 


- 2 v M v M + v v 

” ■ irk R » (E °’ - 


(25) 


with the requirement that < r) < . The inverse of the transformation is now 

applied to obtain from equation (19) 


®$ 0>(X ' E) exp( -°j x) 6 jM 5[ X+R j (e)-r m (e q ) ] 


(26) 


and from equation (22) 


*( 1 ), 

(x,E) =- m, u o, 


1^7 V» vhr exp (- 7 °i t* - V B > - 1 


■ 7 ^[ x + R j( E > + T)* ] 


(27) 


so long as 


^[r m (E 0 ) - x] < R.(E) < ^ R m (E 0 ) - X 


(28) 


Otherwise $P^(x f E) is zero* After a complicated but straightforward manipulation, 
a similar result may be obtained from equation (23) for $j 2 ^(x,E) . 
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The dose as a function of depth is given as 


D(x) 


= 2^ / S - (E) (x,E) dE 

• Jc\ ^ 3 


which is evaluated from equations (26) and (27) as follows: 


( 0 ) 


(x) = A m exp(-o M x) J 6[r m (E)+x-R m (E 0 )] dE 




M 


(tf [ r m <e o ) - x ]} 


exp(-a..x) 

M 


and 


■>'”<*> -E 




a. v. f E u r 
— — I 3 exp /- — 0 . [x - 

V M - "j 1 I 2 5 1 


r^e) - n’ ] 


- <fc[* + V E) + n’ij 


dE 


where 


span the energy limits associated with equation (28) • The 


E and E Q 

11 j 

values of the integrand of equation (31) at the limits of integration are 


exp(-c^x) (E = E^ ) 


exp(-OjX) 


-'v 




(E = V, 


With a linear approximation of the exponent, equation (31) becomes 


D (1) (x) 




A. v. V E JL 

r 3 _ q .) x t eXp( ~ q j x) ‘ eX P (_ °M X) ] 


jM M v M - Vj (a M .j 


( 29 ) 


(30) 


(31 ) 


(32) 


(33) 
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Note that assuming the stopping power to be negligible results in an energy- 
independent solution for the secondary flux: 


-L 


m jM a M 


exp(-a.x) - exp(-a x) 
j M 

a M ' 0 j 


( 34 ) 


which provides some confidence in the above approximations (see also ref. 5). 

In similar fashion, the inverse transformation can be applied to equation (23) 
and D^^(x) can be evaluated as 


D <2> (x, 


A. v. (E' - E' \ 

= 3 3 I u j *jj 

m jk c ^\M a M ( v - v, ) ( a.. - a.)x 

,, M K M 1 

Dk 


exp(-cr.x) - exp(-a M x) exp(-OjX) - exp(-a k x) 


a, - a. 

M ] 


°k - 


(35) 


where 



— R M (E ) - x 
V. M 0 

1 


and 

R j( B L) ■ "5 t[ r h <e o ) - *1 


(36) 


(37) 


The results of equations (36) and (37) are understood to be zero whenever the right- 
hand sides are negative. The above expressions can be applied to various shield 
materials of uniform composition. Each specific application requires knowledge of 
the appropriate transport coefficients Sj(E), c^, and nu^. The next section pre- 
sents the transport coefficients for water. Water is particularly important since it 
resembles organic materials, is under extensive study at heavy ion facilities, and is 
related to an astronaut's self -shielding factors. 
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TRANSPORT COEFFICIENTS OF WATER 


Stopping Power 

In passing through a material, an ion loses the larger fraction of its energy to 
electronic excitation of the material. Although a satisfactory theory of high-energy 
ion-electron interaction is available in the form of Bethe's theory utilizing the 
Born approximation, an equally satisfactory theory for low energies is not available. 
Bethe's high-energy approximation to the energy loss per unit path is given as 


S 

e 


2 4 

4tiNZ 1 Z 2 e 


mv 


In 


2mv 


(1 - p )I. 



(38) 


where Z 1 is the projectile charge, N is the number of targets per unit volume, 

Z 2 is the number of electrons per target, m is the electron mass, v is the pro- 
jectile velocity, p = v/c, c is the velocity of light, C is the velocity- 
dependent shell correction term (ref. 11), and I 2 is the mean excitation energy 
given by 


z 2 ln I 2 ^ f n ln En 


(39) 


where the f n are the electric dipole oscillator strengths of the target and the 
E n the corresponding excitation energies. Note that the sum in equation (39) 
includes discrete and continuum levels. Empirically it has been observed that molec- 
ular stopping power is reasonably approximated by the sum of the corresponding empir- 
ically derived "atomic" stopping powers for which equations (38) and (39) imply 


Z In I 


n .Z . 
3 3 


ln I, 
3 


(40) 


where Z and I pertain to the molecule, Zj and Ij are the corresponding 
atomic values, and nj are the stoichiometric coefficients. This additivity rule 
(eq. (40)) is usually called Bragg 1 s rule (ref. 12). 

Sources of deviations from Bragg's additivity rule for molecules and for the 
condensed phase are discussed by Platzman (ref. 13). Aside from shifts in excitation 
energies and adjustments in line strengths as a result of molecular bonding, new 
terms in the stopping power appear due to coupling between vibrational and rotational 
modes. Additionally, in condensed phase, some discrete transitions are moved into 
the continuum, and collective modes among valence electrons in adjacent atoms produce 
new terms in the absorption spectrum to be dealt with. Platzman proposed that the 
experimentally observed additivity rule may not show that molecular stopping power is 
the sum of atomic processes but rather demonstrate that molecular bond shifts for 
covalent bonded molecules are relatively independent of the molecular combination. 

On the basis of such arguments, Platzman suggested that ionic bonded substances 
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should be studied as a rigid test of the additivity rule because of the radical dif- 
ference in bonding type. He further estimated that ionig bond shifts could change 
the stopping power by as much as 50 percent. Recent results on molecular bond shifts 
on mean excitation energies are discussed in references 14 to 16. Effects of the 
physical state have likewise been studied (ref. 17). 

The electronic stopping power for protons is adequately described by equa- 
tion (38) for energies above 500 keV for which the shell or "tight binding" correc- 
tion C makes an important contribution below 10 MeV (ref. 18). For proton energies 
below 500 keV, it is well-known that charge exchange (electron transfer) reactions 
alter the proton charge over much of its path, so that equation (38) is to be under- 
stood in terms of an average over the proton charge states. It is normal to intro- 
duce into equation (38) an effective charge that is found by averaging equation (38) 
over the charge states so that the effective charge is the root-mean-square ion 
charge and not the average ion charge. At any ion energy, charge equilibrium is 
established very quickly in all materials. Utilizing the effective charge in equa- 
tion (38) appears to make only modest improvement below 500 keV, an indication pre- 
sumably of the failure of this theory based on the Born approximation (ref. 7). 
Customarily data below 500 keV are treated on an empirical basis (refs. 7 and 18). 

The resultant stopping power for protons in water are shown in comparison to the 
evaluated data of Bichsel (ref. 19) in figure 1. 



Figure 1 .- Calculated and experimental stopping powers in water for typical 
cosmic ray ions as a function of kinetic energy. 


The electronic stopping power for alpha particles requires terms in equa- 
tion (38) of higher order in the projectile charge Z^ resulting from corrections to 
the Born approximation. The alpha stopping power cannot be related to the proton 
stopping power through their effective charges. Parametric fits to experimental data 
are given by Ziegler in reference 20 for all elements in both the gaseous and con- 
densed phase. 
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The electronic stopping powers for heavier ions are related to the alpha stop- 
ping power through their corresponding effective charges. The effective charge sug- 
gested by Barkas (ref. 21) is used: 


Z* = z[l - exp(-125p/Z 2//3 ) ] 


(41) 


where Z in equation (41) is atomic number of the ion. 

At sufficiently low energies, the energy lost by an ion in nuclear collision 
becomes important. The nuclear stopping theory used in this paper is a modification 
of the theory of Lindhard, Scharff, and Schiott (ref. 9). The reduced energy is 
given as 


e 



32.53A 1 A 2 E 

+ A 2 H z f 3 + z 


2/3 i 1/2 
2 1 


(42) 


where E is in units of keV/amu and A 1 and A 2 are the atomic masses of the pro- 
jectile and target. The nuclear stopping power in reduced units (ref. 20) is 


r 


1 .593e 


1/2 


(e < o.oi) 


S n- 


1.7 s 1 / 2 ln( e + exp 1) 


1 + 6.8e + 3.4e 


3/2 


(0.01 < s < 10) 


(43) 


ln(0.47e) 

2e 


(10 < e) 


and the conversion factor to units of eV/(10 15 atoms/cm 2 ) is 


f = 


8.426Z 1 Z 2 A 1 


( A 1 + a 2 X Z ? /3 + Z 2 /J ) 


,2/3-v 1/2 


(44) 


The total stopping power Sj is obtained by summing the electronic and nuclear con- 
tributions. Other processes of energy transfer such as Bremsstrahlung and pair pro- 
duction are unimportant. 
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For energies above a few MeV per nucleon, Bethe's equation is adequate provided 
that appropriate corrections to Bragg's rule (refs. 14 to 16), shell corrections 
(refs. 7, 11, and 18), and effective charge are included. Electronic stopping power 
for protons is calculated from the parametric formulas of Andersen and Ziegler (ref. 
18). The calculated stopping power for protons in water is shown in figure 1 in 
comparison with data given by Bichsel (ref. 19). 

Because alpha stopping power is not derivable from the proton stopping power 
formula using the effective charge at low energy, the parametric fits to empirical 
alpha stopping powers given by Ziegler (ref. 20) are used. Applying his results for 
condensed phase water poorly represented the data of references 22 and 23. Consider- 
ing that physical state and molecular binding effects are most important for hydrogen 
(ref. 14), the water stopping power was approximated by using the condensed phase 
parameters for hydrogen and the gas phase parameters for oxygen (which are known 
experimentally). These results are compared with experimental data for condensed 
phase water (refs. 22 and 23) in figure 1. It appears that Ziegler overestimated the 
condensed phase effects for oxygen. 

Electronic stopping powers for ions of charge greater than 2 are related to the 
alpha stopping power through the effective charge given by equation (41). For water, 
the condensed phase formula of Ziegler for alpha particles gives probably the best 
stopping powers for heavier ions. Calculated results for 1 ^0 and ^Fe ions in 
water are shown in figure 1 in comparison with the Northcliffe and Schilling data 
(ref. 24). Good agreement with Northcliffe and Schilling for 56 Fe ions is especially 
important since their data seem to agree with the range experiments of J. H. Chan in 
Lexan 1 (ref. 25). The stopping powers in Lexan and tissue equivalent material can be 
calculated in a way similar to the procedure given above for calculating stopping 
powers in water. 


Nuclear Absorption Cross Section 

The nuclear absorption cross section is calculated from a quantum mechan- 

ical model of the heavy ion reaction. Approximate solution of the coupled-channel 
equations for high-energy composite particle scattering is used to calculate the 
elastic scattered amplitude from which total and absorption cross sections are 
derived (ref. 26). Nucleon-nucleus cross sections are better than 5 percent accurate 
with respect to absorption mean free paths. More accurate nuclear density functions 
are now being incorporated into the model (refs. 27 to 29) and improved calculational 
methods will result in more accurate cross sections at all energies from 20 MeV/amu 
to 30 GeV/amu. The present results derived from reference 26 are shown in figure 2. 


1 


Lexan: 


Trade name of the General Electric Co., Polymers Product Dept. 
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Figure 2.- Calculated (ref. 26) macroscopic absorption cross section 
for typical cosmic ray ions in water. 


Nuclear Fragmentation Parameters 

Nuclear fragmentation by proton beams has been studied by many researchers, and 
approximate formulas for the production cross sections of any fragment with Z > 2 
for proton beams and arbitrary target nuclei have been derived (ref. 30). Basic 
fragmentation studies for carbon ion beams allow extension of the production cross 
sections for Z > 2 and hydrogen beams to any arbitrary nucleus (ref. 30). Given 
the production cross section o(Z^,A^,A^,Z . ,A . ,E) for a fragment of type Aj 
and Zj produced by an ion of type A k afid ^Z k colliding with a target A«p, then 
the fragmentation parameter, or multiplicity, required for the present method is 


V E) ■ 


0( Z k , A^ ,A,p,Zj , A ^ f®) 


W E) 


(45) 


In the present calculation the production cross sections are averaged over the shield 
material constituents. Fragmentation parameters for 50 V at 1 GeV/amu in water are 
shown in figure 3. 
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Fragment charge 

Figure 3*- Calculated fragmentation parameters (multiplicity) of various charge 
states for vanadium fragmentation in water, E = 1 GeV/amu. 


RESULTS 


The dose in MeV/g has been calculated as a function of depth in water for a 20 Ne 
ion beam of 760 MeV/amu according to equations (30), (33), and (35). The calculated 
dose indicates that of the incident 15.2 GeV associated with the kinetic energy of 
each particle, only about half would be ultimately deposited. It is clear in compar- 
ing D 1 (x) with D^ 2 ^(x) in figure 4 that the series approximation to the dose is 

rapidly converging, so that higher order terms are not the source of error. Testing 
for charge and mass conservation of the fragmentation parameters indicates that the 
average total charge of the first reaction products is 4.9e and the average mass is 
9.2 amu. Clearly, considerable mass and charge as well as energy are lost in the 
Silberberg-Tsao fragmentation cross sections. Although the calculational methods 
developed here are converging at least in the case shown, considerable effort must be 
made to improve the description of the fragmentation process. 
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Figure 4.- First three terms of series approximation to dose as a function 
of penetration depth in water for a neon ion beam of 670 MeV/amu. 


CONCLUDING REMARKS 

A suitable approximation of heavy ion transport appears to have been found on 
the basis of the convergence of the series approximation of a 20 Ne beam in water. 
Inaccuracies in the transport coef ficients, particularly the fragmentation parameter 
(multiplicity), do not yet allow meaningful comparison with experiment. Improvement 
in the accuracy of fragmentation parameters would allow eventual calculation of 
important radiation protection quantities with good accuracy. 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
May 12, 1983 
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SYMBOLS 


Aj atomic mass of type j ion, amu 

c speed of light, cm/sec 

C shell correction 

D(x) energy absorbed per unit mass at x f MeV/g 

e electronic charge, coulumbs 

E ion kinetic energy, MeV/amu 

E q incident beam energy, MeV/amu 

E n dipole transition energy, eV 

f n electric dipole oscillator strength 

5 — 1 

Fj ( E) incident beam flux, (cm -sec-MeV/amu) 

I mean excitation energy for stopping, eV 

m electron mass, eV/c 

mj^(E) multiplicity of type j ions produced by collisions of type k ions of 
energy E 

n^ stoichiometric coefficient 

N density of scattering centers, cm - ** 

Rj(E) continuous slowing-down range of type j ion of energy E, cm 
Rj^tRjCE)] inverse function of Rj(E) 

S e electronic stopping power, MeV/cm 

S j ( E) stopping power of linear energy transfer (LET) due to interaction of type 
j ion with orbital electrons of transport medium, MeV/cm 

j (E) specific stopping power or linear energy transfer (LET) due to interaction 
of type j ion with orbital electrons of transport medium, MeV/amu-cm 

S n nuclear stopping power, MeV/cm 

v speed of passing ion, cm/sec 

2 

x one-dimensional position vector, g/cm 

Zj atomic number of type j ion 

* 

Z effective charge of a moving ion 
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p 


ratio of projectile speed to speed of light 


Vj range scale parameter for type j ion 

Cj/'Hj characteristic coordinates of type j ion, cm 

a. (E) macroscopic absorption cross section for type j ion of energy E, cm' 


(x,E) 

differential flux of 
(cm^sec-MeV/amu)" 1 

type 

j 

ions at 

x with energy 

E, 

W 5 : 

) flux of type j 

ions 

in 

Vb 

characteristic 

-3 -1 

space, cm -sec 1 

x j ( 

) flux of type j 

ions 

in 


characteristic 

-3 -1 

space, cm -sec 


Subscript: 

M type of ions in monoenergetic beam 

Superscripts: 

(0 ),(!),( 2) terms in series approximation to equation (13) 
Primes indicate a variable of integration. 
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